Size effects on generation-recombination noise 
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We carry out an analytical theory of generation-recombination noise for a two level resistor 
model which goes beyond those presently available by including the effects of both space charge 
fluctuations and diffusion current. Finite size effects are found responsible for the saturation of the 
low frequency current spectral density at high enough applied voltages. The saturation behaviour 
is controlled essentially by the correlations coming from the long range Coulomb interaction, ft is 
suggested that the saturation of the current fluctuations for high voltage bias constitutes a general 
feature of generation-recombination noise. 



Generation recombination (GR) noise is due to fluctu- 
ations in the number of free carriers inside a device and 
is associated with random transitions of charge carriers 
between states in different energy bands. Accordingly, 
it represents an important noise source in semiconductor 
devices where carrier concentration can vary over many 
orders of magnitude. Typical examples of transitions are 
those between the conduction band and localized levels in 
the energy gap, conduction and valence bands, etc. GR 
noise can be detected as an exoe^_noise when measuring 
current or voltage fluctuations.EJia Lra 

The calculation of the number fluctuations, and hence 
of the current and voltage fluctuations, due to GR pro- 
cesses constitutes in general a rather complicated prob- 
lem. Usually, one has to resort to simplifying assump- 
tions as for instance the local space charge neutral- 
ity approximation and/or the neglect of the diffusion 
current .ErElThe physical implications of these approxima- 
tions remain still an unsolved problem, since an exact 
theory including simultaneously both effects is lacking in 
the current literature. 

The aim of the present paper is to address this issue 
by providing an exact analytical solution of the low fre- 
quency GR noise properties of a two-level resistor model 
within the framework of the drift-diffusion transport the- 
ory. The main result of the exact solution is the evi- 
dence that size effects on GR noise can lead to the satu- 
ration of the associated current spectral density at high 
enough applied bias. Remarkably, present results confirm 
unexpected findings of recent numerical studies by Bo- 
nani and Ghioned concerning uniformly n— doped silicon 
samples with band-to-band and trap assisted generation- 
recombination processes, where the low frequency spec- 
tral density of the voltage fluctuations was shown to 
become independent of the applied bias for moderately 
strong applied bias. 

We consider a two terminal structure consisting of a 
uniformly n— doped semiconductor sample of length L 
terminated by (metallic) ohmic contacts with a single 
trap (donor) level controlling the electron concentration. 



We assume that the transversal dimensions of the sam- 
ple are much greater than the longitudinal dimension, 
thus allowing for a one dimensional electrostatic treat- 
ment. Without loss of generality, only GR noise sources 
are considered. The physical model necessary to describe 
the noise properties of this system consists of the con- 
tinuity equations for the free electrons and the ionized 
traps, the current equation for the electrons, the Pois- 
son equation and appropriate boundary conditions. The 
continuity equations for free carriers and ionized traps 
densities readMcl 

dn(x,t) ldJ n (x,t) 



dt q dx 
dN+(x,t) 
dt 



r n (x,t) +g n (x,t) +7(x,i), (1) 
r n (x,t) +g n (x,t) +j(x,t), (2) 



where n(x, t) and N^(x, t) refers to the free electron and 
ionized trap densities at point x and time t, respectively, 
q is the positive electron charge, J n (x,t) the conduc- 
tion electron current, r n (x, t) the recombination rate and 
g n (x,t) the generation rate. Furthermore, j(x, t) is the 
Langevin noise source related to GR processes-jwhich has 
zero mean and low frequency spectral densitytrau 



-t-oo 

2 / dtj(x,t)j(x,t') 



[9n( x ) +r n (x)]S(x - x'), (3) 



with A being the sample cross-sectional area. The bar 
denotes average with respect to fluctuations. In the drift- 
diffusion approach, the conduction current writes: 



J n (x, t) = q/j,n(x, t)E(x, t) + qD, 



dn(x, t) 
dx 



(4) 



where /i is the electron mobility, E(x, t) the electric 
field and D n the diffusion coefficient satisfying Einstein 
relation, D n — fi n ksT/q, with ks being Boltzmann's 
constant, T the temperature and where non-degenerate 
statistics is assumed. To perform a self-consistent calcu- 
lation, we take into account Poisson equation 



1 



dE{x,t) q 



dx 



= 2. [N+(x,t)-n(x,t)] , 



(5) ±£ 



where e is the semiconductor dielectric permittivity. The 
total electric current is thus given by 



I(t) 



dx 



J„(x, t) + e 



dE{x,t) 

Ft 



(6) 



Finally, as appropriate boundary conditions we take 
ohmic contacts made of ideal metal-semiconductor inter- 
faces in the diffusion approximation, i.e. we assume that 
the diffusion velocity is much smaller than the contact re- 
combination velocity, thus allowing to neglect any contri- 
bution from thermionic emission processes^ Accordingly, 
we take 

n(0) = n{L) ee n c ; Sn(0, t) = Sn(0, t) = 0, (7) 
where the carrier density at the contact is 



(8) 



with Nc being the effective density of states in the con- 
duction band and 4>b n the contact barrier height. Here, 
Sn(x, t) refers to the free electron density fluctuation. 

The model presented above admits an homogeneous 
stationary solution in the form 



n(x) = N t (x) — n c ; E(x) 



qAfin ' 



(9) 



Note that being the density of traps N t and the contact 
density n c two independent variables, one can always find 
for them values which satisfy simultaneously Eq. (||) and 

g n (n, N*) — f„(n, N t ). Under these homogeneous con- 
ditions the system behaves as a resistor with current- 
voltage characteristics satisfying Ohm's law, I = V/R, 
with resistance R = L/(qAfM). 

To calculate the low frequency number fluctuations 



S N (0) = 2 
we use that 



SN(t)SN(t')dt', 



SN{t) = A [ dx5n(x,t), 
Jo 



where 5n(x,t) for low frequencies satisfies 
d 2 Sn 1 ddn 1 



9a; 2 Le dx 



T 2 



-Sn 



T 2 

D 



7- 



(10) 



(11) 



(12) 



Equation (TJ) can be derived by combining Eqs. (|]); 
(^) and (^) after linearization around the homogeneous 
stationary state and by neglecting the time derivatives 
due to the low frequency assumption. The parameters in 
Eq.(n2h are defined as 



Le 

L qV 



L, 



. 9 2 »(l + f)' r n r, 
Here, Ld is a renormalized Debye screening length and 

(14) 
(15) 

From the solution of the second order differential equa- 
tion (JfJ) we calculate SN(t) through Eq. (|TT|), and after 
some lenghtly but straightforward algebra we arrive at 
the following expression for 5at(0), 



f 


_ dr n (n,N+) 


dg n (n,N+) 




dn 


dn 


1 


dr n (n,N+) 


dg n (n,N+) 


n 


0N+ 


dN+ 



where 



and 



S N (0) = 5^(0) - s^io), 



S%(0)=4ALg n r 2 = 4(AN 2 )T, 



(16) 



(17) 



Sf?(0) =4 (AN 2 ) 



(e A ^ 



1) (< 



A 2 L 



1) (Ai - A 2 ) 



2iAiA 2 (A 2 + A 1 )(e^ i - e^ L f 
A a (l - eS Xl+X2)L - 3e A2L + 3e Ali ) (18) 

A 2 (l - e (Al+A2)L - 3e Ali + 



In Equation ( Eq ) Ai and A 2 refer to the eigenvalues of the 
differential operator in Eq. (p"2|) and are given by 



Al,2 — 



1 



2L } 



(19) 



According to standard results of GR noise theory,0 
Si (0) is conveniently written as 



( =) [5j?(0) - 5|?(0)] . 



(20) 



Equation (20), together with Eqs.(|l7|) and (Esh, consti- 
tutes a fully analytical and exact solution for the low 
frequency current spectral density of the problem in sub- 
ject and represents the main result of the present paper. 
In deriving it we have not made any assumption regard- 
ing either the condition for local space charge neutrality 
or the relevance of the diffusion current. We note, that 
in spite of the minus sign in Eq. (|2p|), 5/(0) is a positive 
definite quantity, as should be. 

To investigate the physical properties of the solution, 
we note that the term 5|f (0) vanishes for an infinitely 
long sample thus implying 



5/°(0) = 4 (AN 2 ) ■ 



(21) 



2 



as known from simpler theories.EJfl Therefore, 5|f (0) con- 
stitutes a size effect related to the finite nature of the 
sample. This size effect is analyzed below by investi- 
gating the dependence of the current spectral density 
upon applied voltage. The results are summarized in 
Fig. [I] where we find convenient to normalize the cur- 
rent spectral density 5/(0) to the spectral density cor- 
responding to an infinite sample with thermal current 

It = k B T/{qR), i.e. 5/ T = 4 (AN 2 ) r (l T /N) 2 , and the 
applied voltage V to the thermal voltage Vr — ksT/q. 
With these renormalizing factors, 5/(0) /5/ T (0) as a func- 
tion of qV/(ksT), depends only on the dimensionless 
length, L/L D . 

The essential result we want to stress from Fig. [I] is 
that the spectral density saturates for high applied volt- 
ages to the value 



5|°*(0) = 




The critical voltage for saturation V sa t satisfies 

r y'iO— ; for M- < 1 
Vsat= { , , x2 T L ° , (23) 

I \LnJ 1 ho ~ 

(and an intermediate behavior for f < L/Lo < 10). 
The low bias behavior of 5/(0) displays also interest- 
ing features: it equals 5^°(0) for L/Lo > 10 and 

< V < (l/Lo^J ksT/q, but it is suppressed to 

1/120 (l/Ld^ 5f*(0) for L/L D < 1 and < V < 
V^hO-^p- (and an intermediate behavior in the range 

1 < L/Ld < 10). 

The saturation of 5/(0) and the dependence of the re- 
sults on the ratio L/Lo can be explained under the as- 
sumption that both space charge fluctuations and diffu- 
sion current play a relevant role. Indeed, although the 
system is neutral in average, i.e. n(x) = N t (x), instan- 
taneously space charge fluctuations can appear for which 
5n(x,t) ^ 6Nt~(x,t). The relevance of these fluctuations 
is determined by their interaction with the rest of charges 
in the system through the long range Coulomb interac- 
tion. In this line of reasoning for L < Lo the long range 
Coulomb interaction does not play any role at all, and 
hence, as soon as the applied voltage is higher than the 
thermal value, the low frequency current spectral den- 
sity saturates. On the other hand, when L > Lo the 
long range Coulomb interaction is able to restor space 
charge neutrality for low applied bias. However, as long 
as the applied bias is high enough as to make the tran- 
sit time tt = L 2 /(iiV) shorter than the renormalized 
dielectric relaxation time t c i = e/[g/xn(l + Tt/r n )], the 
long range Coulomb correlations vanish and the low fre- 
quency spectral density saturates.El Furthermore, the fact 



that the diffusion current is playing a relevant role is seen 
on the fact that no shot noise is obtained in the high 
bias regime .as. it is obtained when the diffusion current 
is neglectedocl 

It is worth noting that the results presented in the 
present paper are qualitatively similar to the numerical 
results recently reported by Bonani and GhioneB for a 
more complex GR model. This fact indicates that satu- 
ration at high applied bias can be a general feature of GR 
noise. The present work, then, offers a more complete 
understanding of the finite size effects on generation- 
recombination noise in semiconductor devices. 
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FIG. 1. Spectral density of current fluctuations associated with GR n( 



3 



